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ABSTRACT 

Long-lived high-energy (> lOOMeV) emission, a common feature of most 
Fermi-LAT detected gamma-ray burst, is detected up to ~ 10 2 s in the short 
GRB 090510. We study the origin of this long-lived high-energy emission, us- 
ing broad-band observations including X-ray and optical data. We confirm that 
the late > 100 MeV, X-ray and optical emission can be naturally explained via 
synchrotron emission from an adiabatic forward shock propagating into a homo- 
geneous ambient medium with low number density. The Klein-Nishina effects 
are found to be significant, and effects due to jet spreading and magnetic field 
amplification in the shock appear to be required. Under the constraints from 
the low-energy observations, the adiabatic forward shock synchrotron emission 
is consistent with the later-time (t > 2s) high-energy emission, but falls below 
the early-time (t < 2s) high energy emission. Thus we argue that an extra high 
energy component is needed at early times. A standard reverse shock origin is 
found to be inconsistent with this extra component. Therefore, we attribute the 
early part of the high-energy emission (t < 2s) to the prompt component, and 
the long-lived high energy emission (t > 2s) to the adiabatic forward shock syn- 
chrotron afterglow radiation. This avoids the requirement for an extremely high 
initial Lorentz factor. 
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1. INTRODUCTION 

Gamma-ray bursts (GRBs) are the most luminous explosions in the universe. Their basic 
scenario based on the emission from extremely relativistic outflows with bulk Lorentz factors 
10 2 — 10 3 and isotropic energies of 10 48 — 10 55 erg has been tested, although many questions 
remain open. The Fermi satellite has advanced our knowledge of GRBs significantly, while 
raising some new puzzles. During its first ~ 2yr of operation as of July 27th, 2010, Fermi has 
observed 19 GRBs with photons detected in the LAT (Large Area Telescope) instrument. 
These observations reveal three new properties ( |Granot et al]|2010t |Abdo et al.||2009d| c 



2010| |2009b|a| |Ackermann et al.|2010a|b| |Fermi-LAT fc Fermi-GBM collaborations||20~lT| ) : i) 
A delayed high energy emission, e.g., in GRB 080916C, GRB 081024B, GRB 090510, GRB 
090902B and GRB O9O926A0 ii) A temporally extended high energy emission: at least 10 of 
the first 19 Fermi LAT GRBs have long-lived high energy emission, lasting much longer than 
the burst duration in the sub-MeV band (which declines very rapidly); in 4 out of 10 GRBs, 
the long-lived LAT light curves have a relatively steeper slope, for example, —1.33 ±0.08 for 
GRB 080916C, -1.70±0.08 for GRB 090510, -1.40±0.06 for GRB 090902B, -2.05±0.14 for 



GRB 090926A according to Zhang et al. (2010). iii) A deviation from a pure Band spectral 
function, showing an extra component in GRB 090510, GRB 090902B, GRB 090926A. These 
are among the brightest bursts, but the observations are compatible with the hypothesis of 
having such a component also in the other, less bright bursts, where it is harder to detect 



(Granot et al. 2010). 



Among these 19 Fermi LAT GRBs, GRB 090510 is a short, hard burst, with a duration 



T 90 = 0.30 ± 0.07s( |De Pasquale et aL||2010[ |Ukwatta et aL]|20Q9[ >, located at a redshift 
z = 0.903 ± 0.003 (iMcBreen et al.||2010l |Rau et all 120091). It has been detected by Fermi 



dGuiriec et aT]|2009t [Ackermann et al.||2010b[ ), AGILE ( |Longo et aL"1|2009[ ), Swift ( |Hoversten 



et al. 2009), Konus-Wind (Golenetskii et al. 2009) and Suzaku (Ohmori et al. 2009). Thus 



a large amount of high-quality broadband information is available on this burst, including 
optical, X-ray, MeV and GeV emission. The Swift BAT instrument triggered on GRB 

23 : 00.4 UT, May 10th, 2009 ([Hoversten et al.||2009|, while the 



090510 at T BAT 



00 



GBM instrument onboard Fermi triggered on at T F = 00 : 22 : 59.97 UT, May 10th, 2009 



(Abdo et al. 2009b). Thus, there is a deviation between the two trigger times, which is 
AT = Tq AT — Tq = 0.43s. Hereafter we adopt the BAT trigger time T BAT as a natural 
start time T for computing the afterglow evolution, this being the onset of the main burst. 



The high energy emission of GRB 090510 has all three of the new features we summarized 



1 The first > lOOMeV photons arrive later than the first lower energy photons detected by GBM (Gamma- 
ray Burst Monitor). 
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above: i) the bulk of the photons above 30 MeV arrive 253 ± 34ms later than those below 
lMeV QAbdo et al.||2009b| ); ii) the high energy emission above lOOMeV shows a simple power 
law decay lasting 200s with a temporal decay index «lat = 1-38 ± 0.07 (De Pasquale et al. 
2010[ ) ^ iii) the time-integrated spectrum from T + 0.07s to T + 0.57s is best fit by a Band 



function and a power-law spectrum (Abdo et al. 2009b); the extra power-law component 
photon index of —1.62 ± 0.03 can fit the data well up to the highest-energy (31GeV) photon 



(Abdo et al. 2009b). 



The XRT observations alone give a spectral index 0.57 ± 0.08 (Hoversten et al. 2009), 



while a detailed analysis of the temporal XRT emission combined with the LAT emission 



indicates a spectral index /3 X ranging from 0.51 to 0.81(De Pasquale et al. 2010), and a 
temporal decay index a Xi i = 0.74±0.03 before a break time tx,b = 1.43ks, which subsequently 
steepens to «x,2 = 2.18 ± 0.10. The optical emission initially rises with a temporal index 
« op t,i = — 0-5^oj3, and after a break time t op t,b = l-SS^^ks it decays with a temporal index 

C*opt,2 



1.13±g:io d De Pasquale et al.||2010|). 



Most of the models aimed at explaining the long-lived high energy emission of GRB 
090510 have favored the view that the high energy photons arise from the afterglow emis- 
sion, being generated via synchrotron emission in the external forward shock (e.g., Kumar 



fc Barniol Duran] [20091 [20101 |Corsi et al.||20T0l |De Pasquale et al.|[20T0| |Gao et al.|[2009 



Ghirlanda et al.||2010 


Ghisellini et al. 


2010 


Wang et al. 


2010 


Razzaque 


2010 


) . This expla- 



nation is fairly natural, since an external forward shock model can account for, at least in 
its gross features, not only for the observed delay of the > lOOMeV photons, which corre- 
sponds to the deceleration time-scale of the relativistic ejecta, but also for the long lasting 
> lOOMeV emission, which can be attributed to the power-law decay of the synchrotron 



external forward shock emission (Kumar & Barniol Duran 2009, 2010). However, various of 



the above cited authors use somewhat different readings of the publicly available spectral 
and temporal slope data, which lead them to favor different explanations for the rapid decay 
of the long-lived high energy emission, falling into five different classes of models as follows. 

One set of models interprets the LAT emission as synchrotron emission of electrons 
accelerated in a standard adiabatic ISM forward shock ( Kumar &: Barniol Duran||2009 2010 



Corsi et al. 2010 De Pasquale et al. 2010). These authors argue that they can explain the 
LAT, X-ray, and optical data with plausible parameter values (and we revisit these arguments 
below). However, these authors did not perform a complete enough study to confirm whether 
the whole LAT data including the first 1 second can be explained by this type of models. 
Kumar & Barniol Duran| (|2009l |2010|) fit the late LAT data as F v cx t" 12 , but do not give 



2 Here we use the convention F v oc t a v 
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much significance to the early-time LAT data. Corsi et al. (2010) fit the whole LAT data 
13 ; which seems to explain the early-time LAT data but exceeds the 2a upper limits at 
100 — 1000s). De Pasquale et al. (2010) suggest a steeper decay slope of the 



as f 

late times 

high energy emission (evolving as t~ 1,3S ) for the whole LAT data by taking a larger electron 
distribution index p = 2.5, but they don't take into account in sufficient depth the early-time 
LAT data (t < Is), which is necessary in order to conclude what is the origin of the entire 
LAT emission (see more discussion in Section 6). 



Other models, e.g., Neamus (2010) attribute the high-energy photons to synchrotron- 



self-Compton scattering (SSC) from an adiabatic forward shock propagating into a wind-like 
medium; this, however, requires an extremely small magnetic energy fraction e# = 10~ 10 . 
Another, different adiabatic forward shock model analyses in greater detail the Klein-Nishina 



(KN) effects on the high energy inverse Compton process (Wang et al. 2010). For some rea- 
sonable parameters, the KN effect, as it weakens in time, results in the synchrotron high 
energy emission being increasingly suppressed by the SSC cooling, which steepens the syn- 
chrotron high energy emission decay slope by a factor as large as 0.5. A fourth model views 
the high-energy emission as decaying proportional to which is interpreted as being 

caused by synchrotron emission of electrons accelerated in a forward shock in the radiative, 
rather then adiabatic, regime. For this, the electron population must be significantly en- 
riched, which is attributed to pair production between back-scattering photons and prompt 



outward-going photons (Ghisellini et al. 2010). This model explains the high energy emis- 



sion without considering the constraints from lower energy, e.g., GBM band, XRT band and 
UVOT band emission, and the pair formation becomes inefficient at shock radii larger than 
10 16 cm, while i?dec ^> 10 16 cm for this burst (see §4). A fifth type of model for the afterglow 



of this burst (and others) is a hadronic model (Razzaque 2010), which explains the high 
energy emission as proton synchrotron emission, while attributing the low energy emission 
to electron synchrotron emission from a forward external shock. This requires a large total 
kinetic energy -Ek,iso = 2 x 10 55 er^J with a low radiation efficiency and an extremely small 
fraction of electron energy e e = 10~ 4 . 

In this article, we re-examine the first set of models (the standard adiabatic forward 
shock model) in significantly greater detail than hitherto. We present detailed arguments 
indicating that it is most likely that the forward shock synchrotron emission can only explain 
the LAT emission from ~ 2 — 3 sec. This conclusion disfavors the external shock origin of the 
early-time LAT emission ( Ghisellini et al.||2010 Kumar fc Barniol Duran||2009 ), but supports 



the suggestions that it is related to the prompt emission (Corsi et al. 2010 De Pasquale 



Similarly large energies are required for hadronic models of the prompt emission of this burst (Asano 



et~aTf20l)') 
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. In § 2, we examine the XRT and UVOT observations, and set up a model of the 
long-lived emission based on synchrotron emission from electrons accelerated in a forward 
shock in a uniform ambient environment, including these X-ray and optical/UV observations. 
We discuss the impact of the Klein-Nishina effects on the high energy emission, under the 
constraints imposed by the lower energy observations, which are found to be significant 
for suppressing the SSC cooling. In §3 we use a semi-analytical model to calculate the 
development of the dynamical quantities of the forward shock across the deceleration time 
and into the self-similar phase, and use this to calculate the radiation properties of the long- 
lived high energy emission produced by synchrotron emission. We find a reasonable set of 
parameters which can explain most of the late afterglow, except for the six earliest LAT 
data points in the light curve. In §4 we check several possibilities for the origin of this early 
high energy emission at times t < 2 — 3s. In §5 we discuss the possibility of the line-of-sight 
prompt emission as the origin of the early-time high-energy emission. In §6 we discuss our 
conclusions concerning the most probable origin of the high energy emission from the short 
GRB 090510. 



2. Forward shock model 

2.1. Constraints from low energy emission of GRB 090510 

The afterglow emission of GRB is generally well explained by synchrotron emission from 
electrons accelerated by the shock produced during a spherical relativistic shell colliding with 
an external medium. From the spectral index 0x ~ 0.51 — 0.81 and the light curve slope 



«x,i = 0.74 ± 0.03 (De Pasquale et al. 2010), the closure relation for the X-ray afterglow 



«x,i- 1.5/3x ranges from —0.51 to 0.01, suggesting a slow-cooling ISM external forward shock 
model with < z/x < v\ q which implies that the decay index of the X-ray light curve 



before the jet break is «x,i = 3(p— l)/4 and the X-ray spectral index is f3x = (p— l)/2 (Sari 



et al. |[l998 ;) J^j From the spectral index we can get a constraint on p which is p ~ 2.0 — 2.6. In 



the external shock model, the break seen in X-ray light curve can be explained as a jet break. 



Hereafter we use the subscripts or superscripts 'f and V to represent the quantities of the forward- 
shocked and reverse-shocked regions, respectively, and we use the convention Q x = Q/1Q X in cgs units 
throughout the paper. 

5 The closure relation ax.i — l-5/3x = —0.5 indicates the slow/fast-cooling ISM/wind external forward 
shock model with vx > v\ and the spectral index /3x = p/2, but this implies p — 1.6 < 2.0, which is not 
favored by numerical simulations of shock acceleration (e.g., Achterberg et al. (2001)) or by observational 



data of general GRB afterglows ( Freedman fc Waxman||200T ) . 
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We assume that the jet expands sideways (Sari et al. 1999). In this case the X-ray light curve 
slope steepens gradually from — 3(p— 1)/4 to — p, so that we have 3(p— l)/4 < a x ,2 < P, which 
is reduced to 2.2 < p < 3.9 by taking «x,2 ~ 2.2. Combining the above two constraints we 
have 2.2 < p < 2.6. Then the X-ray decay slope before the break should be 0.9 < ax,i < 1-2, 
which can be acceptable if we take into account the observed fluctuation of the X-ray flux, 
as the light curves we will show below. 

For the rising portion of the optical light curve before the break time, it is natural to 
assume that the optical band is below u^, which induces an optical light curve slope of 
t 1 / 2 (Sari et al. 1998 Kumar & Barniol Duran 2010). The predicted spectral slope z/ 1//3 



is 



consistent with the observations within the large error bars QDe Pasquale et al.||2010[ ). Thus, 
one can try to explain both the X-ray emission and the optical emission before the break time 
with synchrotron emission of electrons accelerated in an adiabatic external forward shock 
with the assumptions z/ opt < u f m and < ux < vl- Then, assuming that the optical band 



is in the regime of z/„ 



m < ^o P t < K 



after the break time, the slope of the optical post-break 



light curve is the same as that of the X-ray light curve. 

In the analytical calculations of this section, we assume that the self-similar phase 
conditions have been established, and for simplicity we neglect the structure of the shock 
wave, considering a spherical shock with a total isotropic energy E and a Lorentz factor T. 
At late times, the adiabatic dynamical evolution of the spherical shock is in the Blandford 
Sz Mckee self-similar phase, where E = I671T 2 R 3 nm p c 2 / '17 is constant and the scaling law of 



the shock wave is T oc R- ?J I 2 (Blandford & McKee 1976b. The shock propagates a distance 



5R ~ 2T 2 c5t/(l + z) during the small observing time St ( Sari||l997 ), and integrating this and 
using the scaling law, one obtains t — (1 + z)R/8T 2 c. 



According to Sari et al. ( 1998 ), the cooling and minimum Lorentz factors of the electrons 



in the forward shock depend on the total isotropic kinetic energy E, the number density of 
the external environment n, and the fraction of the electron energy and magnetic field energy 
e ef and e#f, which can be expressed as 

67r(l + z)m e c 



7 C 



and 



7n 



a T B' 2 Tt 

p — 2m 



2.0 x 10 8 e^. 



-4-^53.5 



3/8 -5/8,1/8 



n 



4 ^2 



[1 + ^)] 



T = 9.0 x 10 4 e ef ,. 



771I/0 — J- J u 1 

-0.4-^53.5^-4 l 2 



1/8,-3/8 



9 P , 



respectively, where 
B' 



1.4 x 10 



-2 1/2 nl/8 3/8 - 

fc Bf,-4- C/ 53.5' i -4 t 2 



3/8 



G 



1 



(2) 



(3) 



(327re#fr nm p c 

and g p = 3(p — 2)/(p — 1), p is the power-law index of the electron energy distribution, y(7c) 
is the Compton parameter of the electrons with Lorentz factor x, and <tt is Thomson cross 
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section. The cooling and minimum frequencies of electrons are 

q e B' 
2Trm P c 



h *t = h^-hlfT = 2.3 x 10 9 [1 + F( 7c f )] 2 e B t^ 2 n-_lt-^eV (4) 



and 



h»l = h^-^llfT = 5.0 x 10 2 e ^_ , e i / f 2 _4^3 / 5t 2 - 3/2 ^eV, (5) 



respectively. The peak flux density of the forward shock synchrotron emission is 

f _ (i + z)jv e vygT _ 1/2 1/2 

^max - 12^~rf| ~ m -4^53.5^-4 ^ J Y> 



(6) 



where the total number of the electrons that the forward shock swept up is N l e = ^nR 3 n 
1.5 x 10 51 £?53 ^n^t^ 4 , and d L is the luminosity distance. 

We get two constraints from the UVOT and XRT data as follows: 



(i) The optical flux density is about 20/xJy at t ~ 100s (De Pasquale et al. 2010), which 
indicates that 

1 /3 

F opt (t = 100s) ~ F^ ax (^P) = 28.5e-^L4t4 E 53 6 5 n^4gp 2/ VJy ~ 20/iJy. (7) 

(ii) The X-ray flux density is about 20/iJy at t ~ 100s, which indicates 

F XRT (t = 100s) ^ Fi ax ^^Xy (P " V2 = 35^i M 6^ 4 ^ 4 n^aTVy - 20/dy. 

(8) 

Combining the above two equations ^ and ^ , we can express the fraction of the magnetic 
field energy and the number density as 

e B{ ,-4 ~ 0.66e~ f *_ 04 £^ 5 #~ 4 ^|, (9) 

n_ 4 ~ o.eee^.o.^^^ 4 ^ 8 - (10) 

where g e i = e p ~ 2 ' 5 and g e2 = e^ p ~ 2 ^^ 3p ~ 1 " > and we can naturally get e ef > O.OAE^^g^gp 1 
from the condition e#f < 1. 

Inserting equations (|9| and (10) into equation ([3]), the downstream magnetic field 
strength is thus constrained to be 

o ~ 10e ef _ 4 ^ 5 3 5 t 2 p p g e2 mG. (11) 

If the shocks involve a magnetic field amplification factor /g > 1 (in addition to shock 
compression), the upstream magnetic field strength would be 

fi u = ^7(417*) ~ itfE&JavG. (12) 
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Thus, the upstream magnetic field strength could be < 5 — 10/zG (e.g., Kumar & Barniol 



Duran 2010), apparently compatible with shock-compression of the typical magnetic field 
in the interstellar medium, with no need of magnetic field amplification. However, as we 
discuss in £j6j the external density deduced here is much below the average interstellar value, 
and likely so would be the external magnetic field, so that additional field amplification may 
be needed. 

Inserting equations (|9| and (10) into equations Q, ^ and ([6]), the characteristic ener- 
gies and the peak flux density of synchrotron emission are therefore 

hv[ ~ 6.6 x 10 9 6^_ . 4 ^ 5 2 3 . 5 t 2 - 1/2 [l + r(7 c f )rW°eV, (13) 

hu\ a ~ 4.0 x 10 2 t~ 3/2 2 e 6 2 °eV, (14) 

and 

FL X ~ 90^/iJy, (15) 

respectively. The above equations show that u opt < until t ~ 1.4ks, and z/^ < u x < v c 
under the constraint e e f -0.4^53.5 > 1.2 x 10^ 3 t^ 4: g~ 1 g^2 , which is easy to satisfy, where 



V(7c) C 1 as discussed in £2.2 Thus, they are consistent with our previous assumptions. 



The radiative efficiency in the assumed slow cooling regime is (Sari & Esin 2001) 



p-2 

Vr = = 6-3 x 10^0^^^^^^, (16) 



which is much less than unity, consistent with our previous assumption of an adiabatic 
forward shock model. 

In order to check whether the LAT emission predicted by synchrotron emission from an 
adiabatic forward shock can explain the LAT observations we need to study two situations. 

(i) Under the condition e e f _o .4^53.5 > O^- 1 ^ , we yield u f m < lGeV < v[ at t < 100s, 
the average synchrotron flux density from the forward shock in the LAT band (lOOMeV to 
4GeV) is 

Hat ~ FL* (^P) ^ ~ 1-4 x 10"%" 3( ^ 1)/ V 1 6 - Vy- (17) 

which is independent of the electron energy e e {E. In addition, the slope of the LAT light 
curve is also constrained, as «lat = 3(p — l)/4, the same as that of the X-ray light curve. 
For p = 2.5, z/ m < lGeV < v\ at t < 100s is satisfied if e e f -0.4-E53.5 > 0.39, then the 
predicted LAT flux is 1.4 x 10~ 3 /iJy at 100 s, consistent with the observational data within 
the error bars at around 100 s. The corresponding slope of the predicted LAT light curve 
is around —1.125, which can explain the late-time data of the LAT observation. If we 
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take a small electron index, for example, p ~ 2.2, we have < lGeV < v\ at t = 100s 



if e e f,-o.4-E l 53.5 > 0.60. Then the predicted LAT flux can be calculated by equation (17), 
approximated as ~ 8.9 x 10 _3 yuJy at 100 s, which is almost one order of magnitude smaller 
than the observed LAT flux. What's more, the slope of the predicted LAT light curve is 
about —0.9, which is too shallow to explain the late-time LAT observation. Thus, in the 
case z/£j < lGeV < v\ at 100 s, we can exclude the p = 2.2 model. 

(ii) If the electron energy satisfies the condition e e f 5 _o.4-E53.5 < §.Ylg~ x g\ 2 , the lower 
end of the LAT band, 100 MeV, is above the frequency u^, and the predicted LAT flux at 
t = 100s can be calculated as 

^ _ p— i p 



Inserting the condition e e f -0.4-E53.5 < 0.12g p x g\ 2 into equation (18), we see that F[ AT < 
4.3 x lO-^g-f-H' (3p ~ 2)/4 /iJy- For p = 2.5, the predicted LAT flux at 100 s is F£ AT (100s) < 
4.3 x 10~ 4 /iJy, which is almost one order of magnitude lower than the observed LAT flux. 
Thus in the case z/lat > v { c at 100 s, we can exclude the p = 2.5 model. 

For p = 2.2, to be consistent with the observed LAT flux at 100 s, the condition 
e e/,-o.4-^53.5 = 6.3 x 10~ 2 is required, constraining the electron energy to a very small value. 
Here we adopt fairly standard values of e e f = 0.16 — 0.6, similar to those observed in long 
GRBs, since these are determined by collisionless shock physics processes on microphysical 
scales, which should be independent of the global properties of GRBs. Then the constrained 
total kinetic energy could be (0.13 — 0.47) x 10 53 erg. Since the isotropic energy at lOkeV — 



30GeV energy band during T + 0.03s-T + 0.53s is E iso = (1.08±0.06) x 10 53 erg (Ackermann 



et al. 2010b), the radiative efficiency of the prompt emission is in the range 70% — 89%. 
Moreover, we constrain e_B/,-4 = SSeJ 3 ^ and n_ 4 = 1.0eg_ 04 by inserting the constraint 



e e/ -0.4-S53.5 = 6.3 x 10 2 into equations ([9j and (10). The constrained i>l is hv l c = 1.1 x 10 7 eV 
at 100s, and v[ < 10 8 eV when t > 1.2s. Since the expected slope of the LAT light curve is 
Q!lat = (3p — 2)/4 = 1.15 for p = 2.2 in the vlat > v i case, similar to that expected from 
the p = 2.5 model in the < z/lat < case, this can also explain the slope of the late-time 
LAT light curve. Even though z/^ decreases by a factor of 1.5 by taking p = 2.2 rather than 
taking p = 2.5, it doesn't change the optical break time by much, and we may not rule out 
the p = 2.2 model in the z/lat > v { c case. 

Thus, both the p = 2.5 model in the z/^ < lGeV < v l c case and the p = 2.2 model in the 
^lat > case can explain the late-time LAT light curve, leading to the same conclusion, 
since the predicted LAT light curves have the similar slope. For presentation purposes, 
hereafter, we discuss the p = 2.5 model in the < lGeV < ul case. 
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2.2. Impact of Klein-Nishina effects on the constraints 



Wang et al. (2010) have studied the Klein-Nishina (KN) effects on high-energy gamma- 



ray emission in the early afterglow, and find that at early times the KN suppression on the 
IC scattering cross section for the electrons that produce the high-energy emission is usually 
strong, and therefore their inverse-Compton losses are small, with a Compton parameter Y 
of less than a few for a wide range of parameter space. This leads to a relatively bright 
synchrotron afterglow emission at high energies at early times. However, as the KN effects 
weaken with time, the inverse-Compton losses increase and the synchrotron high energy 
emission is increasingly suppressed, which leads to a more rapid decaying synchrotron emis- 
sion. This provides a potential mechanism for the steep decay of the high-energy gamma-ray 
emission seen in some Fermi LAT GRBs. 

The Compton parameter for electrons with Lorentz factor 7 e is defined as the ratio of 
the synchrotron self-inverse Compton (SSC) to the synchrotron emissivity, i.e. 



Y(le) = 



SSC 



(7e 



syn(7e 



(19) 



When the KN suppression on the scattering cross section is negligible, Y(j e ) = Y(j c ) is a 



constant for the slow-cooling case (Sari & Esin 2001). However, for high energy electrons 
with a significant KN effect, Y(j e ) is no longer a constant and this affects the electron 
radiative cooling function, as well as the continuity equation of the electron distribution. 
The self-consistent electron distribution is given by 




7m < 7e < 7c 
7c < 7e 



(20) 



for the slow-cooling case (Nakar et al. 2009 Wang et al. 2010), where C\ is a constant. The 



high energy synchrotron photons with energy hv* are produced by electrons with Lorentz 
factor 7* which typically have 7* > max(7 c ,7 m ). Thus, the number density of electrons of 
7* is 



l + ^(7c)^ _ iVs y n(7* 



-Ci7c7* 



(21) 



1 + ^(7*) * 1 + ^(7*)' 

where N syn (%) = Ci 7c [l + Y ( r y c )}% p ~ 1 is the number density of electrons of 7* when only 



the synchrotron cooling is considered (Sari et al. 1998). Therefore, the number density of 



electrons with Lorentz factor 7* is a factor of 1 + Y (7*) lower than that in the case where only 
the synchrotron cooling is considered. Thus the synchrotron luminosity is correspondingly 
reduced by the same factor. We have ^(7*) oc t 1 ^ 2 , in the slow cooling regime, as long as 
7, > Tm e c 2 /hv l m ( |Wang et al.||2010| ). In that case, if y(7*) ^> 1 the synchrotron luminosity is 
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suppressed by the factor V(7*) which is in proportion to t 1//2 , i.e. the light curve decay of the 
high energy synchrotron emission could be steepened by a factor 1/2 at most. Meanwhile, 
the distribution of electrons which are in the region 7 m < 7 < 7 C is not affected by ^(7), 
due to equation 20, thus the lower energy synchrotron emission decay slope is normal. 

The critical photon energy above which the scattering with electrons of energy 7 e just 
enters the KN scattering regime is defined as hiy^(j e ) = rm e c 2 /j e , i.e. 

/m*n(tJ) = rm e c 2 /7c = 0.52 e ; f 2 _ . 4J E; 5 3 1 5 t 2 V + ^(7c)]eV. (22) 

which is much smaller than z/^ under the condition e 2 f _ 4.E53.5 > 0.13t4g~ 2 g 2 2°[l +Y( , yl)]. In 
this case the synchrotron-self Compton scattering is strongly suppressed due to KN effects, 
and 



Y ( 7c f ) [1 + Y ( 7c f 



-4 n 1 w„,fM 7 / 3 4/3 



1.4 x 10" 4 [1 + Y ( 7c f )] 7/3 eZoAK 3 fi /6 9 1 p /3 9l 2 1 - 3 . (23) 



From the above equation, we find that Y^l) <C 1 under the condition e ef _ 04 £?53.5 ^> 
0.023t 5 J 8 gl^g~2' 98 , which is easy to be satisfied before t < 10 4 s. Therefore, since 1^(7*) ^ 
Y(j c ) ^ 1; the distribution of electrons which contribute to the high-energy synchrotron 
photons doesn't change, according to equation (20), which means that the decay slope of 
synchrotron high energy emission can not be affected by Y 



3. Forward Shock Synchrotron Emission Evolution 



We adopt the following parameters for the calculations. From equation (17), we see 



that the flux density in the LAT range is independent of the total isotropic kinetic energy 
and the electron energy fraction at t < 100s, under the constraint e e f -0.4E53.5 > 0.39. Thus, 
we can fix the total energy of the afterglow as E = 3.0 x 10 53 erg, indicating a radiation 
efficiency as 25%. Even if we choose other values of the total kinetic energy under the 
constraint e e f _o.4-£/53.5 > 0.39, our conclusion will not change. Then the fraction of electron 
energy is constrained as e e f > 0.16. We adopt fairly standard values of e e f > 0.16. From the 
equations ^ and (10), we can get the values of e#f = 2.6 x 10~ 3 — 1.3 x 10~ 5 , corresponding 
to an external density range [n = 1.7 x 10 _6 cm~ 3 — 3.3 x 10~ 4 cm~ 3 ]. These densities range 
from much less than the intergalactic medium (IGM), up to a low density inter-cluster 
medium (ICM) or possibly galactic halo baryon density. The critical Lorentz factor which 



is the boundary between the thin and thick shell cases (Kobayashi et al. 2007) is r c 

-,1/8 -1 

y 53.5 n -4 ^-0.5 



[3(1 + z) 3 E/327rnm p c 5 T 3 ] l / s ~ 5.2 x lO 3 ^^:^ 8 ^ 3 / 8 , where the duration of GRB 090510 
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is T = 0.30 ± 0.07s. For an initial Lorentz factor which is not too large, it is reasonable 
to assume that the initial Lorentz factor is smaller than the critical Lorentz factor, which 
means the shell has given the ambient medium an energy comparable to its initial energy 

1 /3 

at the deceleration time tdec = [3(1 + z) 3 E /32irc 5 nm p T ] > T, indicating the thin shell 
case. Then the initial Lorentz factor of the forward shock can be expressed as a function of 
the deceleration time, which is T = 3.3 x lO^^nl^ 8 ^^ 8 . If, for example, we take an 
initial Lorentz factor as suggested by previous Fermi analyses of 2.2 x 10 3 — 5.2 x 10 3 , and 
a not too small number density of n ~ 10~ 4 cm~ 3 , we would obtain a deceleration time of 
^dec ~ 0.3s — 3s. 

To calculate the dynamics of the evolution of the blast wave including the transition 
from the quasi-free expansion, through the deceleration and into the self-similar phase, we 
use the relativistic hydrodynamics equations for the evolution of the shock radius R, the 
mass m swept up by the shock, the opening half-angle of jet 6 and the Lorentz factor of the 
shock T (e.g., Huang et al.||2000"b ) (see Appendix A), and solve these equations numerically. 



The solution of these equations provides the dynamical quantities which we use to calculate 
the radiation spectrum and the light curves discussed in the following. We take into account 
also the jet evolution as it goes through the jet break and starts to expand sideways. The 
half-angle of the jet evolves as dOj = (c s dt'/r) = [(1 + (3)Tc s /r]dt, with a spreading velocity 



in the comoving frame approximated by the sound speed c s (Rhoads 1997, 1999 Huang 



et al. 2000a). After the inverse Lorentz factor becomes larger than the initial jet angle, the 
spreading of the jet speeds up the shock deceleration significantly, which leads to a steeper 
light curve. 

In figures 1 and 2 we present model light curves in the LAT, XRT and UVOT bands 
for two choices of parameters. The forward shock light curves (solid black lines) are shown, 
for each density choice, for two different Lorentz factors and jet opening angles. A smaller 
opening angle is chosen for a larger Lorentz factor in order to get a certain jet break time. 

In figure 1, we show a solution with a density of n — 1 x 10~ 6 cm~ 3 , which would 
correspond to a sub-average density IGM environment. Under the assumptions u opt < 
before the optical break time and v l m < < u ii the steep decay of the latter parts of 
the X-ray and optical light curves can be explained using such a forward shock synchrotron 
emission going through a jet spreading phase, as shown in figure 1 and 2. Note, however, 
that, as in other analyses too, the late-time optical data are challenging to fit. This leads 
to a model forward shock light curve and spectral fit to the X-ray and optical observations 



satisfying all the constraints of £24^ for a set of parameters e e f = 0.17, E — 3 x 10 53 erg, 
p = 2.5, e Bf = 7.0 x 10~ 3 , and 6 = 0.45° for T = 4500, or 6 = 0.38° for r = 9000. The 
collimation corrected energies are 4.6 x 10 48 and 3.3 x 10 48 erg. 
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In figure 2 we show the forward shock solutions for a more moderate density of n = 
10~ 3 cm~ 3 , corresponding to a galactic halo or interarm medium, which satisfy the low energy 
constraints as well as the afterglow epoch LAT, XRT and UVOT data points. The parameters 
are e e f = 0.6, e#f = 1 x 10~ 5 , p = 2.5 and E = 3 x 10 53 erg, for a choice of 8 = 1.3° and 
Tq = 1900, and also for a choice of 8q = 1.1° and To = 3800. The collimation corrected 
energies are 3.9 x 10 49 and 2.8 x 10 49 erg, which is a typical energy for a short GRB. Despite 
the larger e e f in this fit, the forward shock is still in the adiabatic regime according to equation 

In Figures 1 and 2 the LAT energy band light curves predicted by the forward shock 
synchrotron model shown in black solid lines (thicker for the larger Lorentz factor and smaller 
opening angle choice, thinner for the smaller Lorentz factor and larger angle choice). These 
appear to explain well the late time LAT emission, i.e. after 3 seconds, either with the larger 
or smaller angle/Lorentz factor, no matter how small the deceleration time is, but they 
fall below the first six LAT data points. This difference between the early time observed 
LAT data and the synchrotron forward shock emission suggests that there may be another 
radiation component (whose contribution can be represented by the gray dotted lines in 
Figures 1 and 2, corresponding to the high latitude emission of the prompt emission with 
variability timescale At = 0.5s), in addition to the forward synchrotron contribution (shown 
as the black solid lines), at least for its decaying portion. The gray solid lines are the sum 
of the steep decay component and the forward synchrotron emission with a larger Lorentz 
factor (the first, lower data point could be part of either a rising portion or a variable portion 
of the gray dotted component). Thus, it is possible that the early-time high-energy emission 
is not from the afterglow forward shock synchrotron emission, instead having a different 
origin. In the next section we discuss possible origins for the first six data points in the LAT 
band around the deceleration time. 



4. Other Possible Components around the Deceleration Time 

Besides the forward shock synchrotron emission, around the deceleration time there are 
several other possible emission components which could contribute to the early LAT obser- 
vations. Their importance can be estimated using approximate values for the characteristic 
quantities at the deceleration time. At the deceleration time, the Blandford-McKee self- 
similar solution is not yet applicable. For this we can take the usual value of the energy as 



E = ^-R 3 T 2 nm p c 2 and R = 2r 2 ct/(l + z) as the radius at the deceleration time (Sari et al. 



1998), from which we obtain the initial Lorentz factor of the shock and the deceleration 
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radius, which are 



T = 3.3 x lO^^nlf t-f 
R, PC = 3.5 x l^El^nZ'jHlt 



n_ 4 t dec cm. 



(24) 
(25) 



Before we examine the other possible components around the deceleration time, we 
confirm that the forward shock synchrotron emission cannot explain the observation at early 
times in more details. We can then estimate the approximate high energy emission at the 
deceleration time under the previous constraints provided by the low energy observations. 
Here we take the Lorentz factor and radius at the deceleration time as ~ T /y/2 and 
Rd ~ -Rdec/v 7 ^, which are close to the values indicated by the numerical evolution of the 
previous section. Inserting the above Td and Rd into the first expressions of equations Q — 
Q, and meanwhile adopting equations ^ and (10), the characteristic frequencies and peak 
flux density of the forward shock synchrotron emission at the deceleration time are then 



m,d 



6.8 x 10 5 t- e f^ 2 °eV, 



ao 2 



hvl d ~ 3.9 x lO-e^cM 



7712 / 2 „2 -1(1 v 

^53.5*dcc 9p9e2 eV > 



F , ~ 

max,d 



i2^ 2 2 VJy- 



(26) 

(27) 
(28) 



Ghisellini et al. (2010) suggests the radiative forward shock model to explain the steep 
temporal decay of high energy emission, which requires > v\ to get the fast cooling case. 
However, the radiative model cannot explain the shallower decay at later time (t > 10s) (as 



seen in Figure 7 in Ghisellini et al. (2010)), which agrees better with an adiabatic model in 



slow cooling case at that time. According to the constrained characteristic frequencies, the 
fast cooling case at the deceleration time requires that 



e ef < 1.7 x lO-^t^Vs" 



(29) 



which is inconsistent with their assumption of a very high energy fraction of electrons e e f 



0.9. This is because Ghisellini et al. (2010) assume the late-time emission to be also in the 



fast cooling case, and do not take into account the constraints from the low energy emission. 



In addition, using equations ^ and (10), we get an estimate for the forward shock 
synchrotron emission at lGeV at the deceleration time, 



lGeV,d 



max,d U<J 



p-1 

2 



0.049t 



-3(p-l)/4„-2.6 
dec 



9ei ' ^Jy- 



(30) 



This is almost one order of magnititude below the observed LAT flux density, for the case 
when the deceleration time is smaller than 3s (this is seen also in Figures 1 and 2). The 



- 15 - 



synchrotron forward shock emission cannot explain the early-time high-energy emission (the 
first six LAT data points near the deceleration time), although it can explain very well the 
late LAT results (beyond the sixth LAT data point, when the afterglow can be considered 
as established), which are consistent with the numerical results in ^3] 

Thus, we find that there has to be another component, contributing to the early high 
energy emission. In the rest of this section we consider several such possibilities, such 
as a reverse shock synchrotron component, a reverse shock SSC component or a cross IC 
component (forward/reverse shock synchrotron photons scattered by electrons from reverse- 
shocked/forward-shocked region), a high latitude (curvature) component of the prompt high 
energy emission. 



4.1. Reverse Shock Synchrotron Emission 

We consider the thin shell case as is assumed in §3, the flux peaks at the crossing time of 
the reverse shock t x = tdec > T with Lorentz factor T x = Td, and the ratio of the comoving 
number densities of the forward-shocked to that of the reverse-shocked regions is given by 



n T /n { ~ T d ~ 2.3 x lO 3 ^^ t d ^ (Kobayashi et al. 2007). The internal energy densities 



e and the bulk Lorentz factors T of the two regions are equal with each other (Zhang et al. 



2003). Consequently, we have that 

M' 2 Mi §=^ 3 [l + n7e, x )]" 2 , # £ = ?Kb, (31) 



c 

1 



characterized by the ratios = f|^J 2 > 



€er 



The minimum and cooling frequencies of the reverse shock synchrotron emission at the 
crossing time under the low energy constraints are 

/< x = l.le e f,^^ 4 ti /4 »B,i^^ eV, (32) 

K,x = 3.9 x lO^A^S^flJ^ 10 [1 + Y(jI !X )] ~ 2 eV. (33) 

According to §2.2, we have Y(jl) < Y(Y C ) ^ 1 around the deceleration time, thus the 
synchrotron spectrum does not change and the ratio of the reverse shock synchrotron flux 
to forward shock synchrotron flux can be calculated as 



FlGeV.x = F max,x ( \ " 2 ( h K,x \ 2 f nA ^ »p-l«x(p+l)/2 / h K,x \ 

^icev.x UGevj \nj e B {lGeV J 



(34) 
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Inserting equations ^ and (10) into equation (34), the ratio turns to be 

■^IGeV.x _ n 9 „ p/2 p/4+3/2 3p/8-W/2-W-l 6.7 »/2 -p-3.0 /„kn 

f — u.zot ef _ 4^53 .5 t dec ji b1 Ji £t0 y e i g p g e2 ■ {oo) 

r lGeV,x 

The ratio could be as high as 0.59 since the total kinetic energy can be as high as E ~ 
9.0 x 10 53 erg by considering a reasonable radiation efficiency > 10%. Considering equation 



(30), the flux density from reverse shock synchrotron emission cannot be as high as ~ /zJy 
to explain the observations at early times. Thus, the reverse shock synchrotron emission is 
unlikely to contribute to the early-time LAT observation. 



4.2. The SSC and EIC Emission 

Besides the reverse synchrotron emission, we consider the other four possible IC pro- 
cesses, including the synchrotron self-Compton (SSC) processes in forward and reverse 
shocks, and two combined-IC processes (i.e. scattering of reverse-shock synchrotron pho- 
tons on electrons accelerated in forward shocks and froward-shock synchrotron photons on 
electrons accelerated in reverse shocks). 

The optical depth of inverse Compton scattering in the revers e shock in the Th omson 
regime is r r = \a^R d /T d = 7.7 x 10- 8 e^_ ^^ 4 t^ /8 ^ /2 ^ 5 -°dWang et al.||200l|). The 



optical depth to inverse Compton scattering in the forward shock in the Thomson regime 
is r f = |a"Tni?d = 3.8 x l0^ 12 el { _ 0A E^/ 5 2 t]^gpg~2' - Taking the peak flux of synchrotron 
emission at the crossing time of the reverse shock and forward shock, which are -F max = 
OSOe'^l^El^t^^BAgp #e2° J y and f L = 12 9e2°^Y, respectively, we can get the peak 



flux density of the four Inverse- Compton components as follows (Wang et al. 2001) 



= ^ max = 2.2 x 10" 8 Jye 2 ef> _ A t^ B ,i9p9a 2 - , (36) 
^ = -^Lx = 8.8 x 10-^yeV% A E^%f g^g 7 ™, (37) 
*S = - f ^ax = LI x lO^Jyc^A 1 ^^/ 2 ^ (38) 

and 

= ^Lx = 4.5 x 10~ 17 Jy^f.-o^m^te^S", (39) 
respectively, where the superscripts 'rr' and 'ff' mean the SSC emission in reverse shock 
and forward shock respectively, 'fr' and 'rf mean the scattering of reverse shock photons 
on the electrons in forward shocks and forward shock photons on the electrons in reverse 
shocks. Therefore, the IC contributions (except for the SSC emission in the reverse shock) 
can be excluded, since even their peak fluxes are much smaller than the early-time LAT 
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observations. Although the peak flux of the SSC emission in the reverse shock is close to the 
observed flux, the reverse shock SSC emission can also be excluded because its flux peaks 
at a low energy of about 20keV, indicating that the flux in the LAT band is about 6 orders 
lower than the peak flux. A low flux for these IC processes is mainly due to a very low 
circum-burst density in GRB 090510 inferred from the low-energy observations. 



4.3. The High Latitude Prompt Emission 



The spectrum of GRB 090510 shows a power-law component in the LAT band (> 
lOOMeV), whose physical origin at early times is unclear. One might consider the early part of 
the extended emission (shown by the gray dotted line with a decay slope «Lat(0.37 — 3s) = 2.0 



in Figures 1 and 2) is due to the high-latitude emission of the prompt emission (Kumar & 



Panaitescu 2000). Because photons from high latitude regions with respect to the line of 



sight will arrive later than that from low latitude region due to the curved front surface of 
the jet, one observes a fast decreasing emission rather than an abrupt stop of the emission, 
the so-called "curvature effect" . Then the high latitude emission flux of the prompt emission 
evolves as 



F u (t) oc 



At 



(40) 



according to Toma et al. (2009) (The details are shown in Appendix B), where to is the pulse 



peak time. According to Ackermann et al. (2010b) and Abdo et al. (2009b), the variability 



timescale derived from the BGO emission (lOOkeV— few MeV) is about 14 ± 2ms in the 
time interval 0.17 — 0.37s since the BAT trigger, but the variability timescale for the LAT 
emission is not determined. If we assume that the variability timescale for the high energy ( 
larger than lOOMeV) emission, which is dominated by the power law component, is the same 
as that of BGO observation, then we adopt At ~ 0.01s in equation (40). Thus, the high 
latitude flux of the prompt emission at t + 0.1s is F(t Q + 0.1s) = ll _2_/3 F(to)- Since the 



index of the power law component in the time interval 0.17 — 0.37s is 1.66 ±0.04 (Ackermann 
et al.||2010b| |Abdo et al. | |2009b| ), we take the rough value (5 ~ 0.66. Then the flux of the 
high-latitude emission at to + 0.1s will be reduced by a factor of ~ 600 relative to that of 
the line-of-sight emission at to ~ 0.37, which is too steep to explain the early high energy 
emission at t = 0.37 — 3s. 

Since the origin of the LAT emission in the prompt phase, dominated by the extra 
power-law component, is unclear and can be different from that of the Band component in 
the GBM energy range, the variability timescale in the LAT range can generally be different 
from At ~ 0.01s. In particular, the emission radius of the extra high energy component may 
be much larger than that of the Band component, which may lead to a much larger pulse 
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width in the LAT energy range than in the GBM energy range (see e.g., Toma et al. (2010)). 
If we take At ~ 0.5s, then the high latitude flux of the prompt emission at t ~ 0.37s 
evolves as F(t) ~ (2t + 0.26s) _2 ' 66 F(0.37s), which together with the synchrotron forward 
shock emission with a deceleration time tdec < (2 — 3)s can explain the early high energy 
emission at t ~ 0.37 — 3s, shown in Fig 1. 



5. The Line-of-Sight Prompt Emission as the Origin of the Early-Time 

Afterglow High-Energy Emission 

As a final possibility, we consider the possible influence of the tail-end of the prompt 
emission in the description of the early afterglow. In §3 and §4 we assumed a small GRB 
duration T 90 = 0.30 ± 0.07s, following De Pasquale et al. (2010) based on the observations 



in the GBM and BAT bands. Under this assumption, we discussed the forward shock 
synchrotron emission, the forward reverse shock synchrotron emission, the four possible 
crossed inverse Compton processes, and the high latitude emission of the prompt emission, 
and we concluded that none of them can explain the early decaying part of the high energy 
emission at 0.3 < t < 3s except for the high latitude portion of the prompt emission with a 
larger variability time scale plus the synchrotron forward shock emission with a deceleration 
time tdec < (2 — 3)s. 



We focus on the BAT detections around T + (1 — 3) s, shown in De Pasquale et al. 



(2010), which implies that the prompt emission could last until such times. We calculate the 
flux of the synchrotron forward shock emission at t = 3s for p = 2.5 in the BAT energy range, 
and obtain ~ 52yuJy (the flux is similar for p = 2.2). This is about one order of magnitude 
smaller than the observed flux. Thus these BAT detections are not from the external shock, 
but may be attributed to the prompt emission. 



The LAT prompt emission may also last until such times. |Ackermann et aL (2010b) 
obtain no significant temporal correlation between Nal (8keV — 260keV) data and the LAT 
data > lOOMeV between t = 0.17 s and t = 0.47 s. This does not necessarily indicate 
that the LAT prompt emission (the extra component) and the GBM prompt emission (Band 
component) have fully unrelated origins. Some theoretical models suggest that the two 



components may arise at different radii (e.g., Toma et al. (2010), Wang et al. (2006)) in 
the same shells, or the two components may arise by the leptonic and hadronic processes in 
the same internal shock (e.g., Asano et al. (2009); Razzaque ( 2010[ )). Therefore it appears 
reasonable to consider that GRB 090510 has a duration as long as 2 — 3 s and the high 
energy emission before t ~ 2 — 3s is a part of the prompt emission. 



-19- 



Under this new assumption, with a longer GRB duration T ~ 2.0s, one infers a critical 
Lorentz factor T c ~ 2.0 x 10 3 ^{ 8 5 n:J /8 T 7 3 3/8 . For the thin shell case with an initial Lorentz 
factor To smaller than T c , the shell will be decelerated efficiently by the reverse shock after the 
GRB main duration. The high energy light curve with the initial Lorentz factor To = 1900 
and the number density n = 10~ 3 cm~ 3 is shown by the thin solid line in Figure 2, which can 
explain the late-time high energy emission well. If we assume that the LAT emission around 
t = 2s is a flare of duration At = Is, which peaks at t = 2s with a peak flux F(t Q ), and 
having ( De Pasquale et al.]|2010 ) a spectral index (3 = 0.41+qJ® at t — 1.5 — 2.5s, then the flux 



due to curvature effect at t — 3.5s is F(t) = O.HF(to), which is compatible with the seventh 
data point. In this case, the deceleration time could in principle be even larger, implying a 
smaller Lorentz factor, but due to the lack of precise information about the last flare, e.g., 
its duration and the flux peak time, we cannot use this to derive a firmer deceleration time. 

For a larger initial Lorentz factor Tq > T c , the reverse shock will transition from 
the Newtonian phase to the rel ativistic phase at the time £n = {E /lQitAniripC^T^) 1 ^ 

o.o8ii4 /2 5 r -j. 7 n:fs 



Sari 



1997), with the width of the shell in the comoving frame A = 
cT/(l + z) = 3.1 x 10 10 T .3cm. The deceleration time when the dissipated energy is compa- 
rable to the total kinetic energy of the shell is about tdec = 2T. At the time < t < ta cc , 
the Lorentz factor of the forward shock evolves as T oc i -1 / 4 , and the radius of the shell 
evolves as R oc t 1 ^ 2 . Thus, the characteristic frequencies and the peak flux density evolve 
as z/^ in oc t~ l , v x c oc t _1 and F max oc t. As a result, the flux of high energy emission in the 
region u f m < z/ L at < v f c evolves as -Flat oc t<~ p + 3 )/ 2 ; and the flux density of high energy 
emission in the region i^lat > v l c evolves as -Flat oc £(~ p + 2 )/ 2 . This shows that the forward 
shock synchrotron emission in the thick shell case cannot explain the early-time high energy 
emission due to its flatter light curve. And the reverse shock emission in the thick shell case 
can also be excluded, using similar calculations as in §4. 

Based on the above analysis, while the forward shock synchrotron afterglow dominates 
the LAT emission at late times, the early times LAT emission must be attributed to a prompt 
emission, which can be limited from above by a high-latitude curvature emission envelope. 



6. Discussion and Conclusions 

We have investigated whether the photons received by the LAT can be explained solely 
by the forward shock afterglow emission or whether the prompt emission necessarily con- 
tributes to the early part of this emission. We have addressed this question with the help 
of the broad-band observations from the XRT and UVOT instruments onboard of Swift and 
the LAT onboard of Fermi for the case of the short GRB 090510, obtaining constraints on 
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the extended high-energy emission based on the XRT and UVOT observations. 



We have obtained a good fit to the XRT and UVOT observations (except for the late- 
time optical data) in terms of synchrotron emission from an adiabatic forward shock in 
spreading jet model, which constrains the environment of this short burst, implying a low 
number density n ~ 10~ 3 — 10 _6 cm -3 , consistent with a binary progenitor scenario (e.g., 



Belczynski et al. 2006), corresponding to an electron energy fraction e e f ~ 0.6 — 0.2 and a 
weak magnetic field with an energy fraction est ~ 10~ 5 — 10~ 2 . The latter could in principle 
be due to an upstream magnetic field of ~ just shock-compressed, without need for 



magnetic field amplification (Kumar & Barniol Duran 2010). However, the magnetic field 
values of ~ fiG derived in our galaxy correspond to an average ISM of density n ~ 1 
cm -3 , whereas the fits for GRB 090510 indicate external densities comparable to a halo or 
intergalactic medium with n ~ 10 -3 — 10 _6 cm -3 , where based on flux-freezing the field would 
be expected to be much less than nG. Thus, the magnetic field may still need to be amplified 
in the shock. We note also that while the electron energy fraction can be as high as 0.6, the 
adiabatic condition is still satisfied, since the low energy emission constrains the radiation 
to be in the slow cooling case, indicating a radiation efficiency much smaller than unity, in 



contrast to the radiative regime assumption of Ghisellini et al. (2010), where they did not 
take into account low energy constraints. 



In De Pasquale et al. (2010), it is assumed that the LAT data at t — 1 — 200s with 
slope —1.38 arises from synchrotron emission in the regime z^lat > ^ i- e -' v c < 10 8 eV, in 
contrast to our conclusion that the synchrotron emission in the regime v c > lGeV evolving 
with slope ~ —1.1. Part of the difference may be due to De Pasquale et al. (2010) using the 



formulae of Granot & Sari (2002), while we used the formulae of Sari et al. (1998). However, 
applying the constraints on e#,n, e e iE in |De Pasquale et al. (2010) and taking p = 2.5 in 
the expression for z/ from Granot & Sari (2002), one obtains v ~ 6.3 x 10 7 eV at 100 s. 



Therefore, z/ c f = 6.3 x 10 8 eV at Is since v c oc t 1 / 2 , falling in the LAT band, inconsistent 



with the De Pasquale et al. (2010) assumption. The simple power-law decay shown in their 



Figure 1 is not applicable due to the v c crossing through the LAT band, and it overestimates 
the LAT photon flux by a factor of 1.7 at t — Is and 2.1 at t ~ 0.4s . We note also that the 



formulae in Granot & Sari (2002) is applicable only in the self-similar phase, which takes a 



few e-folding times after the deceleration time to fully develop. Our numerical calculation 
shows that around the deceleration time, the flux should be reduced by a factor of ~ 2 — 3. 
Furthermore, if the deceleration time is not larger than the duration T, the LAT light curve 
will be flatter at t < T, which cannot explain the early-time LAT emission, as discussed in 



}5 and also strongly argued by Maxham et al. (2011) 



The initial jet angle is found to be constrained to be 9 S < 1°, which implies a collimation- 
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corrected kinetic energy in the range of short GRBs. We have used for our numerical calcula- 



tions the semi-analytical jet spreading model of Huang et al. (2000b) (detailed in Appendix 
A), assuming that the observer line of sight is approximately along the jet axis. This leads 
us to identify a feature in the late XRT and UVOT light curves as symptomatic of a jet 



break. We note however that van Eerten et al. (2010) argues that numerically the light curve 



break related to a jet edge as seen by an off-axis observer may appear hidden for a longer 
time (weeks) than it would be for an on-axis observer, and the alternative possibility cannot 
be ruled out, although this being one of the brightest bursts reduces the probability of its 
being observed from an off-axis direction. 

We have considered the origin of the high energy emission under the constraints provided 
by the low energy observations. Calculation of the Klein-Nishina effect on the IC scattering of 
synchrotron photons by electrons accelerated in the forward shock indicate that the Klein- 
Nishina effect suppresses the high energy IC emission significantly, which results in a Y- 
parameter much smaller than unity. The distribution of the high energy electrons is not 
affected by the Klein-Nishina effect, so that the synchrotron high energy emission cannot 
steepen due to a suppression of the IC scattering through Klein-Nishina effects as argued in 



Wang et al. (2010). 



The duration of the prompt emission of this GRB is a crucial parameter for establishing 
the physics of the LAT emission. Based on the phenomenon that the Fermi GBM emission 



turns over sharply at t ~ 0.30 ± 0.07s (De Pasquale et al. 2010), we adopted a duration of 



T = 0.30 ± 0.07s, compatible with a thin shell case where the deceleration time is larger than 
the GRB duration. Here, however, we find that a forward shock synchrotron emission model 
agrees well only with the late-time (t > 3s) high-energy observations of the LAT, but this 
model falls well below the early-time high-energy emission (the first six LAT data points). It 
cannot contribute to the power-law high energy component in the prompt spectrum, which 
suggests that the early-time high-energy emission may have a different origin from the late- 
time high-energy emission. We excluded various other possibilities for the early-time high- 
energy emission, such as synchrotron emission from the reverse external shock, SSC emission 
from reverse/forward external shock, and IC scattering of forward/reverse synchrotron pho- 
tons by electrons accelerated in reverse/forward shocks, as well as the high-latitude prompt 
emission with a small variability timescale. The latter is too steep to contribute to the 
early-time high-energy emission, due to the assumed high temporal variability of the prompt 
emission, unless the early-time high-energy emission variability timescale is as large as 0.5 s. 

We are led to the conclusion that the early-time high-energy emission is likely due to 
the final portions of the prompt component, or else to the high latitude component of the 
prompt emission with a variability time scale as large as 0.5 s plus the synchrotron emission 
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with a deceleration time t^ec < (2 — 3)s. This is supported by the fact that we cannot 
explain the first six LAT data points by synchrotron forward shock emission if we consider 
the combined data in the three bands (XRT, UVOT and LAT). In this case, the long-lived 
high-energy emission can be naturally explained as the result of the prompt emission at early 
times (t < (2 — 3)s), and the afterglow emission at late times (t > (2 — 3)s) from the adiabatic 
forward shock synchrotron radiation, with reasonable parameters. In this two-component 
model of the long-lived high energy emission, the shock deceleration time of tdec = 2s results 
in a fit with an initial Lorentz factor around T ~ 2000, not much larger than the lower limits 
Tmin = 950 ±40 and 1220 ±60 (Abdo et al. 2009b) obtained from the pair-production limits 
implied by the presence of 3.4GeV and 31GeV photons. With smaller values of tdec, very 
large initial Lorentz factors are needed (but see Ioka (2010)). Although there is so far no 
statistically significant evidence for an early steeper decay, due to the sparse nature of the 
LAT data, it is striking that by eye there are 4 out 10 GRB long-lived LAT light curves which 
have relatively steeper slopes than the slopes, ~ — 1, of typical late-time X-ray and optical 
afterglows, according to Table 3 and the light curve fits in the figures of Zhang et al. (2010). 
The steep decay slope of the LAT light curve is difficult to explain with a normal external 
shock model. A similar conclusion was subsequently suggested by Liu & Wang (2010), and 
also by Maxham et al. (2011 ). Our two-component model also explains the steep high-energy 
temporal decay index values olat ~ 1.5 as being the natural result of the superposition of 
the tail-end of the prompt regime and the start of the afterglow. While our model for the 
high energy afterglow is in principle applicable to other LAT bursts as well, a homogeneous 
sample of similarly detailed data including low energy constraints will be needed, on a large 
number of objects, before definite conclusions can be reached. 
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A. Jet Deceleration and Spreading Dynamical Model 



We solve the following equations for the evolution of the shock radius R, the mass m 
swept up by the shock, the opening half-angle of jet 9 and the Lorentz factor of shock V 



e.g., Huang et al. 2000b) 



The evolution of the shock radius is described by 

_^cr ( r + Vf-i), 

with t as the observer's time and f3 = \/T 2 — 1 /T. The swept mass evolves as 

dm 



dR 



2ttR (1 — cos6)nm p , 



(Al) 



(A2) 



with m p as the proton mass and n as the number density of surrounding medium. The 
evolution of the opening angle considering the jet spreading is described as 



de 

~di 



1 da 
R~~dt 



3 (T + y/T^T) 
R 



(A3) 



where a is the comoving lateral radius of the jet (Rhoads 1999 Moderski et al. 2000) and 



the comoving sound speed c s = 7(7 — l)(r — 1) 1+ ^ r _^ c 2 with the adiabatic index 7 as 
7 m (4T + l)/(3r). The conservation of the total energy can be expressed as 

dr r 2 - 1 

dm M e j + r] r m + 2(1 — e)Tm ' 



(A4) 



where the radiative efficiency is rj r = for the adiabatic case, and M e j is the ejecta mass. 



B. The High-latitude emission 

Because photons from high latitude region with respect to the line of sight will arrive 
later than that from low latitude region due to the curved front surface of the jet, we will 
observe a fast decreasing emission instead of an abrupt cutoff of the emission, which is so 
called "curvature effect". 



We neglect the radial structure of the emitting shell for simplicity. Under this assump- 
tion we calculate the flux from the shell at radius rj which expands toward us with Lorentz 



factor Ti by (Granot et al. 1999 Woods & Loeb 1999 Ioka & Nakamura 2001 Yamazaki 



et al. 


2003; 


Dermer 


2004 


Toma et al. 


2009 



F v {t) 



1 + z 
di 



[i + r 2 2 (t)p 



(Bi) 



-24- 



where the photon frequency in comoving frame is 



v' = (1 + z)u- 



2I\ 



(B2) 



with v as the photon frequency in observing frame, and 9(t) describes the emitting point of 
the high-latitude emission. We have the comoving emissivity 



^oci/^ocOL + r]^)) 



(B3) 



where /3 is the spectral index of t 



spectra fit of (Abdo et al. 



re observed emission, which is consistent with the power-law 



2009b 



with tj as the emission time in t 



By using the relationship 6(t) = 
re central engine frame, we get 

t - (t - At) 



1/2 



At 



(B4) 



i + r^ 2 (t) 

where t = (l + z)(ii— ^ ) is the observed peak time of the high energy pulse and At = n ^2^ 
is the observed dynamical time scale. Therefore, the observed flux evolves as time as following 



F v {t) oc 



't- (t - At)' 
At 



-2-/3 



(B5) 
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Log10[t/s] 



Fig. 1. — GRB 090510 broad-band fitting. The data are the energy flux densities averaged in 
the observed energy bands (De Pasquale et al. 2010): LAT(100MeV-4GeV, filled squares); 
XRT (0.2 keV - 10 keV, crosses); UVOT (white band, open circles). The black dotted, 
dashed and solid lines are the light curves of the optical, X-ray and high energy emission, 
respectively, from the adiabatic forward shock synchrotron emission in the jet spreading 
model with parameters E = 3 x 10 53 ergs, e e = 0.17, = 7 x 10~ 3 , n = 1 x 10~ 6 cm~ 3 
and electron index p = 2.5. The thick lines are for an initial Lorentz factor T = 9000 and 
an initial jet half-angle 9q — 0.38°; the thin lines are for T Q = 4500 and 6 = 0.45°. The 
thick gray dotted line shows a hypothetical component, which can be explained by the high 
latitude emission of the prompt emission with a variability timescale At = 0.5s, which makes 
up for the difference between the synchrotron emission with T = 9000 and the observations. 
The thick gray solid line is the sum of the thick gray dotted line and the thick black solid 
line. 
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Fig. 2.— GRB 090510 broad-band fitting via the adiabatic forward shock synchrotron 
emission in the jet spreading model, with parameters E = 3 x 10 53 ergs, e e = 0.6, e# = 10~ 5 , 
n = 10 _3 cm~ 3 and p = 2.5. The thick lines are for an initial Lorentz factor To = 3800 and 
an initial jet half-angle 9 = 1.3°, while the thin lines are for T = 1900 and 9 = 1.1°. The 
thick gray dotted line shows a hypothetical component, which can be explained by the high 
latitude emission of the prompt emission with a variability timescale At = 0.5s, which makes 
up for the difference between the synchrotron emission with To = 3800 and the observations. 
The thick gray solid line is the sum of the thick gray dotted line and the thick black solid 
line. 



